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Abstract The binding neuron model (Vidv bidin[l998h is 
inspired by n umerical simulat ion of Hodgkin-Huxley-type 
point neuron dVidvbidalll996h. as well as by the leaky in- 
tegrate-and-fire model dSegundo at all. Ill968l) . In the bind- 
ing neuron, the trace of an input is remembered for a fixed 
period of time after which it disappears completely. This is 
in the contrast with the above two models, where the post- 
synaptic potentials decay exponentially and can be forgot- 
ten only after triggering. The finiteness of memory in the 
binding neuron allows one to construct fast rec urrent net- 
works for comp uter modeling dVidvbidal l2003). Recently, 
(Vidybida 2007), the finiteness is utilized for exact math- 
ematical description of the output stochastic process if the 
binding neuron is driven with the Poisson input stream. In 
this paper, the simplest possible networking is considered 
for binding neuron. Namely, it is expected that every out- 
put spike of single neuron is immediately fed into its input. 
For this construction, externally fed with Poisson stream, the 
output stream is characterized in terms of interspike interval 
probability density distribution if the neuron has threshold 2. 
For higher thresholds, the distribution is calculated numer- 
ically. The distributions are compared with those found for 
binding neuron without feedback, and for leaky integrator. 
It is concluded that the feedback presence can radically alter 
spiking statistics. 
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1 Introduction 

The main function of a neuron is to receive signals and to 
send them out. In real neurons, this function is realized fhrou£ 
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concrete biophysical mechanism, the main parts of which 
are ion channels in excitable membrane and variations of 
ionic concentra tions inside and outside of nerve cell and its 
processes, see (ISchmidtlll975h for details. The same func- 
tion might be realized by means of any other mechanism 
able to support signals processing in the manner, which is 
characteristic of a real neuron. If so, then it would be in- 
teresting to develop a model, which realizes in an abstract 
form a concept of signal processing in real neurons, and is 
exempted from necessity to follow any biophysical mecha- 
nism supporting the processing. Such a model is necessary 
for quantitative mathematical formulation of what is going 
during signals/inform ation processing in neural systems, see 
(vanHemmen 2007) for discussion. Attempts to develop 
such a model are mainly concentrated around concepts of 
coi ncidence detector a nd temporal integrator, see discussion 
in dKonig et al. II1996I) . O ne more model, t he binding neu- 
ron (BN), is proposed in d Vidybida I fl998). This model is 
inspired by numerical simulation of Hodgki n-Huxley-type 
neuron stimulated from many synaptic inputs d Vidybida 1 199i 
as we ll as by the leaky integrate-and-fire model ( Seg undo at all 
1968). It describes functioning of a neuron in terms of events, 
which are input and output spikes, and de gree of temporal 
coher ence between the input events, see (Vid ybida 1 1 19981 
l2007h for details. 

It is observed, that during processing of sensory signals, 
the spiking statistics of individual neurons changes substan- 
tially when the sign al travels from per iphery to more cen- 
tral areas (see, e.g. (lEggermont 1 119911) ). The changing of 
spiking statistics could underlie the inf ormation condensa- 
tion, which happens during perception (Konig and Krugerl 
2006). This transformation of statistics may happen due to 
feedforward and feedback connections between neurons in- 
volved in the processing. Having in mind such possibilities, 
it would be interesting to check what happens with spike 
h train statistical properties when it passes neuronal structures 
with feedback connections. 

Usually, feedback/recurrent connections are considered 
between several neurons. In this paper we consider the sim- 
plest possibility, namely, the single neuron with feedback. 
Such a configuration, which we regard as the simplest possi- 
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ble n etworking, can be found in real biological objects (see, 
e.g. dAroniadou-Anderjaska et al.ll999l : lNicoll and Jahr Il982h ). 
As neuronal model we use binding neuron as it allows to 
obtain exact mathematical expressions suitable for further 
analysis. It is expected that input stream in any synapse of 
the neuron is Poisson one. In this case, from mathematical 
point of view, all inputs can be replaced with a single one 
with Poisson stream in it, having its intensity equal to the 
sum of all intensities in the synapses (Fig[TJ top). The bind- 
ing neuron works as follows. Any input impulse is stored in 
the neuron during time T and then it is forgotten. When the 
number of stored impulses, E, becomes equal to, or larger 
then the threshold one, iVo, the neuron sends an output spike, 
clears its internal memory and is ready to receive impulses 
from the input stream. One obtains the binding neuron with 
feedback (BNF) by immediate feeding each output impulse 
to the neuron's input (Fig[TJ bottom). In this case, just after 
firing, the neuron has one impulse in its internal memory, 
and this impulse has time to live equal T. 

The specifics of mathematical analysis of BN-type sys- 
tems is due to presence in those systems both deterministic 
and stochastic dynamics. Namely, the neuron obtains its in- 
put from a random stream (stochastic component) and every 
impulse is stored for the same fixed period of time (deter- 
ministic comp onent). This is i n the contrast with the mass 
service theory dKhinchinll955h . where the service time (coun- 
terpart of time to live, t) is random, Poisson-distributed. The 
simultaneous presence of deterministic and random dynam- 
ics in real neurons is due to the fact that in real neurons 
the impulse existence in a neuron (exposed as the excita- 
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tory postsy naptic potential) is suppor ted by electrochemical 
transient (Hod gkin and Huxievlll952h . which is determinis- 
tic, whereas the input impulses come from other neurons and 
external media in irregular (random) manneiQ. 



2 Condensing of Information with Neurons 

It is widely accepted that during flow of sensory signals 
in a hierarchical manner from sensory periphery to central 
brain areas, the information, which is present in the sig- 
nals, becomes less analogue and more discrete, eventually 
resu lting in representing disc rete symbols or entities (see 
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t — memory 

Fig. 1 Schematic representation of binding neuron (top) and BN with 
feedback (bottom). 



After processing with the BN, the output spike train con- 
sists of two output impulses at moments and fg- The BNF 
gives three spikes at moments tj, t4 and ?6- It is clear that 
the output trains contain less information then the input one. 
Nevertheless, having the output train for either BN, or BNF, 
one can make some condensed conclusions concerning the 
input train. 

In the case of BN, the output spike at moment f 3 tells us 
that there where input impulses at moments f 3 and ?*, where 
t* e]?3 — V,t-i[, and the input impulse at moment t* did not 
trigger an output one at the moment of its arrival. For real- 
ization shown in Fig. 12 f* = t%. Information delivered in the 
output spike t% is indeed condensed, because the presence 
and exact timing of this output spike does not change if po- 
sition of f * deviates remaining within interval ]t% — T; tj [. The 
same is about output spike at moment t^. 

In the case of BNF, the output spike at moment f 3 tells us 
that there where input impulses at moments / 3 and t*, where 
t* e]? 3 - T;f 3 [, and the time interval ]/* - v,t*[ is free of in- 
put impulses. The same is about t^. Two output spikes at 
moments t$, ?4, where ?4 — ? 3 < T, tell us additionally that in 
the input there were 3 consecutive impulses separated by in- 
tervals shorter then T. Similar conclusions can be made for 
binding neurons with No > 2. Moreover, condensed conclu- 
sion about input, which is based on corresponding output 
spike train, can be formulated for other neuronal models, 
like Hodgkin and Huxley, or leaky integrator. The difference 
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of information within the flow must decrease in order to map 
various input spike trains from the sensory periphery into 
the same discrete entity. This process of consecutive reduc- 
tion of information is known as condensation. We now put 
a question: What could be the primary element in which the 
condensation takes place? It seems that single neuron is a 
suitable candidate for such an element. In the case of bind- 
ing neuron this can be explained as follows. 

Consider an input spike train like upper train in Fig|2j 
The train can be regarded as signal from a receptor neuron. 



tion in precise and clear mathematical manner, whereas for 
other models it does not. A less precise, fuzzy formulation, 
which is suitable for any model is that the output spike sig- 
nals about presence in the input train temporally coherent 
(distributed over short time interval) sets of impulses (see 
( Vidv bida 1119961 Il998l) for discussion). 



3 Output Intensity of BNF with Threshold 2 



1 Compare with lGoel et al.l i Ti^jBritvina and Eggermoml {HI), The intensity of output stochastic process can be defined by 
where cases consistent with mass service theory are analyzed. three different ways: 
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Fig. 2 Example of input spike train and corresponding output for 
binding neuron (BN), and binding neuron with feedback (BNF). In 
both cases No = 2. 



1. It is possible to define the instantaneous intensity (see 
dKhinchin Ill955l) ). X a (t), as the probability to obtain an 
output impulse at moment t in infinitesimal interval 5 di- 
vided by s: 



X [t) — lim • 



(1) 



where w(s,t) denotes the probability to obtain impulse 
in the interval [t;t + s[. We do not intend to analyze the 
X (t) time dependence here. 
2. As limit, or mean output intensity one can choose the 
following 



X a = lim X (t) . 



(2) 



3. In the definition of X in (0, the role of time limit is to 
ensure that the initial state of the system is forgotten. In 
this case one can define/calculate X a as the factor in the 
expression X dt, which gives the probability to obtain an 
output impulse in the infinitesimal interval dt, if nothing 
is known about previous states of the neuron. 

It can be shown that both definitions 2 and 3 bring about the 
same value for X a . Therefore, we choose the third one here. 
Calculations bas ed on the sec ond definition can be fulfilled 
with the help of dFeller II 19661 Part XI, §8). 

The probability to obtain an output impulse from the 
BNF with threshold 2 in the interval dt , X a dt, can be calcu- 
lated as product of probabilities of two independent events: 
(i) an input impulse is present in the dt; (ii) the interval be- 
tween that impulse and its predecessor is not longer than 
X. If the input stream is Poissonian, then the probability of 
event (i) is A dt, and of event (ii) is 1 — e' x T , where A is the 
intensity of input stream. Thus, 



)A. 



(3) 



4 Distribution of Output Intervals for BNF with No = 2 

Let us consider a BNF with threshold Nq = 2 and internal 
memory X, which obtains its input from Poisson stream with 
intensity A. Thus, the neuron fires every time when input 
impulse comes T, or less units of time after its predecessor. 



The output statistics can be described in terms of the 
probability density distribution to obtain an output interspike 
interval t with precision dt. For this purpose it is enough to 
calculate the probability, Ph(t,x)dt, of the following event: 
the next firing happens f units of time later than the previous 
one. Let the input impulses, which come after the previous 
firing, are numbered with numbers 1,2, 

The above-mentioned event can be decomposed into sev- 
eral alternatives, which are numbered with the number k of 
input impulse, which triggers the next firing. Notice, that for 
t < x only one alternative is possible. It happens if the first 
input impulse comes not later then X units of time after the 
previous firing. In this case, the neuron still keeps impulse 
from the previous firing, and the input causes the threshold 
achieving and firing. There is no other ways to get output 
interval t < X. Thus, for t 6 [0;t[, the probability density 
distribution is as follows 



P b (t,T)dt =e Xt Xdt . 



(4) 



It is impossible to obtain output interval t > X with a 
single input impulseQ Thus, for t > X, possible alternatives 
are numbered with numbers 2, 3, . . . , k max , where k max = 
[t/x] + 1 , and [x] denotes the integral part of x. 

Assume, the k-th alternative is realized by input arrival 
times t\ , ?2, • • ■ , fjfc-i, fjfc = t. Not all arrival times are admitted 
(see (O and further, below). In accordance with the defini- 
tion of Poisson process, the probability of such realization is 
given by the following expression: 

e~ Xtl Xdt ie -^ t2 - tl) Xdt 2 ■ ■■e- X{ '-' k - l) Xdt . 

The probability Pbk(t, x)dt that the £-th alternative is realized 
with any admissible values of t\, t2, fjt-i can be calcu- 
lated by integrating of the above expression over the region 
of (k— 1) -dimensional space with coordinates t\, f2, ■ • ■ , ffc-i, 
defined by the following conditions: 



fl > T, fl +T< t 2 , 



_2 + T<f*_i <t. 



(5) 



and t — < X. The required integral over the region de- 
fined by d5j can be calculated exactly: 



t-(k-2)t t-(k-3)x t 
e~ Xt X k - 1 j dt x J dt 2 --- J dt k _ l Xdt = 

.^■- ''-{'-■f' M. (6, 

{k-iy. 

If k = kmax, then (|5]l ensures: (k — l)-th impulse is in the 
interval ]t — X;t[, and &-th impulse at moment t will cause 
firing. Thus, in this case 



P bk (t,x)dt = e 



-Xtik-\ 



-i ('-(*-!)*) 



k-\ 



(*-!)! 



-Xdt , k = k„ 



The value of Pb(t, t) can be chosen arbitrary. 
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Fig. 3 Interspike intervals (ISI) distribution P/,(r, t) for % = 10 ms, X 
= 0.01 ms"', No = 2. Left — calculated in accordance with (|4j, 0, 
right — calculated numerically. 



If k < k max , then integral (O includes also configurations for 
which f£_] < f — t. For these configurations fc-th input im- 
pulse at moment ? will not cause firing. The contribution of 
these configurations into the integral (O is given by the fol- 
lowing expression 



e- Xt l k - x 



r— (*— l)T t-(k-2)x t-X 

j dt\ J dt%--- J dtk_\Xdt 



X <i+T tk-2+X 

6 X (*-!)! Ut < 
which should be subtracted from ©. Thus, for 2 < k < k m 



P bk (t,x)dt = 

X k ~ l 



-Xi 



([t-{k-\)x) k - x -{t-kx) k - l )xdt. 



(*-l) 

The total probability is calculated by summation over all 
alternatives. Notice, that k max changes by 1 when f passes 
through integer multiple of T. Thus, for m = 1,2, . . . the fol- 
lowing statement is valid: if mx < t < (m+ l)x, then 



Q m 

P b {t,x)dt = e- Xt — (t-mx) m ldt- 
ml 



(7) 



-Xt 



2<k<m 



^({t-{k-\)x) k -'-{t-kx) k - l )Xdt. 



For t e [0;t[ the function P b (t,t)dt is given by 0. The 
distribution Pb(t, x)dt is analogous to distribution P(t)dt — 
e~ X dt known for Poisson process. The graph of Pb{t, x) is 
shown in Fig|3] 



5.1 Connection with BN distribution 

It is interesting that function P(t, x), which gives the prob- 
ability den sity distribution f or binding neuron without feed- 
back (see d Vidvbida 1 120071 Eq. (6))) has simple intercon- 
nection with P b (t,x). In order to find this interconnection, 
denote restriction of P/,(t, x) onto interval [mx; (m+ 1)t[ as 
Pbm(t,T). then means: 

mx<t < (m+l)x =*> P b {t,x)=P bm (t,x) 7 m= 1,2,... . 

Substitute here t + x instead of t: 

mx<t + x< (m+ 1)t =3- P b (t + x,x) = P bm {t + x,x), 

where m = 1,2, . . . , or 

(m-l)x <t <mx P b (t + x 1 x)=P bm (t + x,x), 

where m = 1,2,... . Substitute here m instead of (in — 1 ) : 

mx <t < (m+l)x => P b (t + x,x) =P b;in+ i{t + x,x) 1 

in = 0, 1,2, ... . 

The explicit expression for P bM+ i (t + T, t) can be obtained 
from ©: 



Pb,m+iit + t,x) =e 



-X(t+r) 



x m+1 

(m+1) 



■(t-mx) m+l X+ 



\k-l 



I e~^^((t-(k-2)x) 



k-l 



2<k<m+\ 



{t-{k-\)x) k ~ l )X 



The last expression coincides with the corresponding term 
in the (IVidvbidall2007l Eq. (6)) multiplied by e~ Xr . Thus, 
the following representation takes place: 



0<t <x P,,(t,x)=e- x, X, 
x<t ^ P b {t,x)=e- XT P(t-x,x). 



(8) 



The last express ion together with the fact that P(t, x) from 
dVidvbidall2007l . Eq (6)) is normalized, allows one to check 
easily that P b (t,x) is normalized as well: 



/ 



P b (t,x)dt = 1. 



5 Properties of the distribution 

Notice that after firing, the neuron starts from standard state: 
it keeps a single impulse with time to live equal x. There- 
fore, there is no correlation between consecutive interspike 
intervals. 



5.2 Mean interspike interval 

Having for P b (t, x) representation one can easily calcu- 
late mean interspike interval, W, which is defined as 




o 
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Substitute here representation 

T oo 

W-f.rHt + f.^'Pff-rt*- 
} -Te- lT + e-^ J{t + x)P{t,x) 



l-e 



dt 



l-e 



-Xx 



■e~ Xx I tP{t,t)dt. 



The last integral is calculated in dVidvbida 1120071 Sec. 3.2). 
Use found there expression: 



W 



l-e 



-Xx 



+ e" AT - (2 



A \ e^-l 



A(l 



-Xx\ 



(9) 



5.3 Coefficient of variation 



Coefficient of variation for obtained distribution Q can 
be calculated as follows 



W2 
W 2 



(10) 



where W is given in (0, and W2 is the second moment of 
distribution (0: 

w 2 = / t 2 p b (t,z) dt. 


Here use representation ([8]): 

T 00 

W 2 = A J t 2 e~ Xt dt + e~ u J t 2 P(t -z,z)dt= (11) 



2-((Ar) 2 + 2AT + 2)e 
X 2 



-Xx 



+ e~ Xr J (t + x) 2 P(t,T) dt. 



The second term here can be split into three: 

e -Xx J z 2 p ^ t) ^ dt = e- u x 2 , 




e - AT 2T y tP(t,T)dt =e~ AT 2r- I 2+ At ' - 


(used same expression as for calculating (01), and 



-At ' 



1 



J t 2 P{t,T), 



e~ Xz J rPu.x)dt. 
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Fig. 4 ISI distribution P^(t, t) found numerically for T = 10 ms, A = 
0.05 ms~ 1 . Left — No = 4, right — No = 6. Notice the discontinuity at 
t = t. 30000000 triggerings were taken in both cases. 



The ISI distribution's second moment for BN without feed- 
back can be calculated similarly as it is done for its first mo- 
ment. This gives 



J t 2 P(t,x) dt=j^- 



2 le 2Xx + (AT-3)e AT +1 



X T 



(e x *-l) 2 


Substitute this into ( fTTT) . This gives 

2( ,XX e AT + ^ T 

W2 = 



(12) 



(13) 



X 2 {e^-\f' 
Substitute this and (0 into dlOl) . this gives 

c fev = ^Are-^ + l. 
Coefficient of variation gets its maximum value, c bvm , 

Cbvm = V^e' 1 + 1 sa 1.32 

at At = 1 (Fig. |5). 

It is also possible, by using Eq. (1121) . to calculate coeffi- 
cient of variation, c v , for BN without feedback: 



c v = 



/ 2At < ? ;u + 0.5 1 

4 e 2AT_ 4e AT + 1 + 2' 



The c v gets its maximum value equal to 1 at At = 0, and 
decreases monotonically when At increases (Fig. [5j. 



6 Numerical Simulations 

Numerical simulations were executed here for several pur- 
poses. The first purpose was to check numerically correct- 
ness of the expressions found analytically in previous sec- 
tions. A C ++ program was developed, which allows to cal- 
culate the Pi,(t,T). The Poisson streams were generated by 

transformation of uniformly distribu ted sequences of ran- 

dom n umbers (see, e.g. Eq. (12.14) in ( [Computational Cell Biology] 

Those sequences were produced with the system pseu- 
dorandom number generator from libc library in the Linux 
operating system, as well as with th e Mersenne Twister, mt 19937 
dMatsumoto and Nishimura1 fl998). The two methods give 
indistinguishable results. The program includes the BNF class, 
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which analyzes the input stream and fires in accordance with 
the rules, described above. With the help of that class, output 
stream samples were produced by calculating N = 30000000 
output spikes. The samples are scanned for interspike inter- 
vals of various duration, and the probability density distri- 
bution is then calculated by normalization. The numerically 
obtained ISI distributions are in good agreement with analyt- 
ical expression, as it can be seen in Fig. [3] Also, the second 
moment of Pj,(t,x) was calculated numerically for several 
values of A, T. Deviation of numerically found values from 
analytical expression ( 113b is within 0.01% ~ 0.1% range. 

The second purpose of numerical calculations is to ob- 
tain ISI distributions for higher thresholds. The above math- 
ematical reasoning for BNF with threshold No =2 becomes 
extremely cumbersome for thresholds No = 3, or higher. It 
seems that there should be a more effective mathematical 
approach for higher thresholds. Meanwhile, it is possible to 
calculate numerically the probability density distribution for 
any threshold value. Examples of calculated densities are 
shown in Fig. [4] These densities are in qualitative agreement 
with what is found analytically for No = 2, except of the 
fact that the initial part of ISI distribution is increasing for 
No > 2, whereas for No — 2 it is decreasing. The initial (for 
t < X) part of the probability density distribution Pj,(t,x)dt 
can be easily found analytically for any threshold No. In- 
deed, denote the moment of the previous firing as 0. At this 
moment BNF stores one impulse with time to live X. The 
next firing happens at moment t < x iff No — 2 input im- 
pulses come within the interval ]0;?[, and one more impulse 
within [t;t + dt[. The probability of such event for Poisson 
process is known, which gives for any No>2 

Pb{t^)dt = e- X, y >_ Xdt, t<x. 

This function is decreasing for No =2 and increasing for 
higher No, which explains seeming qualitative disagreement 
between No = 2 and No >2 cases. 

The third purpose was to compare the ISI distributions 
found here for the binding neuron model with those for leaky 
integrate and fire (LIF) model. In the program developed, 
the BNF class was replaced with LIF class, which repro- 
duces the simplest version of the LIF model. Namely, the 
LIF neuron is characterized by a threshold, C, and every in- 
put impulse advances by yo the LIF membrane voltage, V. 
Between input impulses, V decays exponentially with time 
constant Xm- The LIF neuron fires when V becomes greater 
or equal C, and V = just after firing. Examples of the ISI 
distribution obtained for various parameter values are shown 
in Fig.|6] 



7 Conclusions 

We calculated here the intensity and output interspike inter- 
vals distribution for binding neuron with feedback, which 
is stimulated with Poisson stream. For BNF with threshold 
Nq = 2 this is done analytically, for higher thresholds — 



numerically. It is interesting to compare the obtained dis- 
tributions with those known for other models. In Fig. [5] the 
distribution is shown for binding neuron without feedback. 
Curve 4 in Fig. [5] is q ualitatively similar to distribution ob- 
tained numerically in ( Segun do at all. I[l968h for leaky inte- 
grator model in a slightly different stimulation paradigm. By 
comparing these distributions with those found here for BNF 
(see also comparison of coefficients of variation in sec. l5.3b , 
one could conclude that even the simplest possible network- 
ing is able to change radically statistical properties of spik- 
ing process. This gives a hint about what could take place 
with spiking statistics of individual neurons in a network. 



1.4 | 1 1 1 1 1 0.4 




0.7 1 1 ~ r 1 ' 1 ' ' 1 1 ' ' —3 

02468 10 0123456 



Fig. 5 Left: Coefficient of variation as function of x = Xz for BN 
(lower curve) and BNF (upper curve). Right: Interspike intervals dis- 
tribution P{t,z) for BN with No = 2 and without feedback (from 
(Vidybida 2007)). Here X = 1 s _I ; curves number 1,2,3,4 correspond 
to x = 0.5 s, 1 s, 2 s, 3 s. 

Numerical calculations made for the LIF model (Fig. [6]l 
suggest, that introducing feedback might result in qualitative 
changing of spiking statistics for other neuronal models as 
well. 



8 Discussion 

The model of binding neuron used here is simplified in a 
sense, that it does not follow time course of ionic currents, 
or transmembran e voltage. The purpose of this model (see 
dVidvbidalll998h ) is to formulate in abstract form the an- 
swer to the question: What does neuron do with signals it re- 
ceives? The question well can b e answered in the framewor k 
of more detailed models, like dHodgkin and Huxlevl[l952h . 
But usage of more detailed models for description of less 
basic functions, like neural coding, or information process- 
ing, would be the same as to describe computer functioning 
in terms of Kirchhoff 's laws: it is correct, but not productive. 

The exact discontinuities in the output ISI distributions, 
which can be seen in Figs. |3|4| are due to abrupt loss of feed- 
back input influence x units of time after triggering. Output 
ISI, which is shorter then x, is created with the feedback 
spike involved. The longer ISIs are created without feed- 
backed spike involvement. Therefore, the jump is in the di- 
rection of smaller probabilities. In the models, in which the 
influence of input spike diminishes gradually, one could ex- 
pect the decreasing region of probability density function 
course in the range, where role of feedback inputs becomes 
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Fig. 6 ISI distribution I%(t) found numerically for leaky integrate and 
fire model with feedback. Used 30000000 output spikes. Firing thresh- 
old, C = 20 mV, input intensity, X = 0.1 ms~'. Left: membrane time 
constant, Tm = 3 ms, input impulse amplitude, yo = 15 mV. Right: 
%M = 6 ms , y Q = 7.5 mV. 



pers are incorporated in a complicated network. Thus, their 
spiking statistics is influenced by other neurons. Therefore, 
a more developed network model is needed in order to com- 
pare with experimental data. Such a model will be studied in 
future. 
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small. This could cause a bimodal distribution of output ISIs, 
like shown in Fig. [(J right. Nevertheless, for special parame- 
ter values, the genuine discontinuity can be as well observed 
for the LIF model, like in Fig.il leffl. 

The coefficients of variation dependence on x = At can 
be explained as follows. For x — > both BN, and BNF out- 
put streams become Poissonian. Consider the BN case. The 
BN will generate an output spike in interval [t;t + dt[ if 
three conditions are satisfied: (i) there is input spike in [t;t + 
dt[, (ii) the previous input was received at t — T, or later, 
(iii) the previous input did not triggered BN. Violation of 
cond. (iii) with (i), (ii) satisfied is improbable when At — > 
0, because this means appearance of two consecutive input 
ISIs, both shorter then T. For Poisson input this may hap- 
pen with probability (l — e~^ T ) , and for small x may be 
neglected. In this case the desired probability of output is 
(l — A dt, which describes Poisson stream with inten- 
sity A' = (l — e T ) A. For this stream, coefficient of vari- 
ation will be 1 . Similar reasoning are valid for BNF. In the 
opposite case, when At — > °°, violation of condition (iii) for 
BN cannot be ignored. Actually, for high stimulation rates, 
the BN will act as perfect integrator. The output stream of 
perfect integrator is y-distributed, with c v < 1 . For BNF at 
high stimulation rates, every feedbacked spike will combine 
with next input one, and trigger next output spikfl This pos- 
sibili ty was mentioned as "dancing in step" in dMacKav I 
1962, p.43). In such a regime, output stream exactly repro- 
duces the input one, hence, is Poisson stream with q, v = 1 . 
For intermediate values of At the "dancing in step" will be 
interrupted from time to time by waiting longer then T for 
the next input spike. The triggering, which is next to this 
event, must happen without feedback involvement. Combi- 
nation of this two possibilities gives maximum variability of 
output stream at A t = 1 . 

Finally, it would be interesting to compare ISI distribu- 
tions found here with those observed experimentally. The 
configur ations with f eedback are known for real biological 
objec ts, dAroniadou-Anderjaska et al J 1999tlNicoll and Jahr I 
The self-excitating neurons described in the cited pa- 

3 For the LIF model, presence of discontinuity in ISI distribution in 
Fig. [6] left, can be proved mathematically rigorously. 

4 iVo = 2 is expected. 
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